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Abstract 

The associated production of a J/0 and a photon in photon-hadron 
collisions is considered, and shown to be a good probe of the pres- 
ence of colour-octet mediated channels in quarkonia production, as 
predicted by the factorization approach within NRQCD. Total and 
differential cross sections for photoproduction at fixed target experi- 
ments and at HERA are presented. Associated J/tp + 7 production at 
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1 Introduction 



The study of heavy quarkonia production is a good testing ground for our 
understanding of the transition region between the realms of perturbative and 
non-perturbative Quantum Chromodynamics (QCD). Different production 
models have been devised, and all of them have to deal with the problem of 
properly describing and matching these two phases of QCD. Most recently 
the so-called Non Relativistic QCD (NRQCD) factorization approach (FA) 
has been proposed by Bodwin, Braaten and Lepage (BBL) an effective- 
field-theory framework in NRQCD is used to separate the short-distance scale 
of annihilation and production of heavy quarkonium from the long-distance 
scales associated to the quarkonium structure. A brief review of the various 
models is presented in section 

In particular the BBL approach (also referred to as FA) has been applied 
to describe quarkonia production at the Tevatron providing a successful 
picture || of the production rates of J/ip and ip' - much larger than earlier 
predictions - which is related to the relevance of colour-octet contributions in 
the production mechanism. It is to be noted that this successful description 
relies on a certain number of non-perturbative parameters having to be fit- 
ted to the data. Therefore additional and more extensive comparisons with 
experimental data coming also from other kind of reactions are necessary to 
finally assess the validity of this approach. To this aim, calculations of in- 
clusive quarkonia production in e + e~ annihilation [|J, Z decays ||, hadronic 
collisions at fixed-target experiments || , B decays j7| and 7p collisions (both 
for S'-wave [Bl-CT and P-wave [0,0] states) have recently been performed 
within its frame. From refs. [^,10 1, in particular, the predictions based on the 
leading colour-octet contribution appear not to agree with recent experimen- 



tal data obtained at HERA by the HI Collaboration | l3fl , indicating either 
a reduced phenomenological importance of the colour-octet terms than sug- 
gested by the Tevatron data analysis || or the possible relevance of higher 
order corrections. 

In this paper we will consider the associated production of a J/0 and a 
photon, where the photon has to be directly produced in the short- distance 
interaction, i.e. it should not come, for instance, from the radiative decay 
of a P-wave quarkonium into a J/ip. This process therefore must have a 
very distinctive experimental signature, namely a J/ip balanced in by the 
hard photon. We will argue that the study of this reaction in 7p collisions 
provides a powerful tool for assessing the importance of colour-octet terms 
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in quarkonia production. First results of our analysis have already been 
presented in jl3| . 

The paper is organized as follows: in section ^| we will briefly review 
various quarkonia production models and the BBL approach. Then in sec- 
tion § we will introduce the basic diagrams which contribute to J/tp + 7 
photoproduction and discuss why this process provides a clear signature of 
colour-octet terms. Numerical results and Conclusions will follow in sects. 4 
and 5 respectively. 



2 Production Mechanisms 

Any approach attempting to describe the production (or decay) of a heavy 
quarkonium must deal with two issues which - for sufficiently large quark 
masses - can be kept distinct: the production (or the annihilation) of the 
heavy quark- ant iquark pair (QQ) constituting the quarkonium and their 
binding into a single physical long-lived particle. The heavy quark mass, 
rriq ^> Aq CD , sets the scale for the short-distance interaction of producing or 
annihilating a QQ pair, which can therefore be studied within perturbative 
QCD (pQCD) fll5| . The binding of the two quarks into a bound state is on 
the other hand a long-distance process, and non-perturbative effects have to 
be included at this point. These are usually parametrized via form factors 
describing the probability of the QQ pair to form the bound state. The de- 
gree of rigour, completeness and complexity of this sector varies greatly from 
approach to approach: 

The Colour Evaporation Model (CEM) rests on duality arguments 



17fl in assuming that cc pairs produced with an invariant mass below that of 
a DD mesons threshold, (taking charm as an example of heavy quark) will 
eventually hadronize into a charmonium state. Being this assumption only 
qualitative this model of course is unable to predict the production rate of a 
particular quarkonium state. It is therefore not very suitable for the study 
of exclusive final states. 

The Colour Singlet Model (CSM) ]TJ] tries to overcome this difficulty 



of the CEM by making a very precise request: the QQ pair must be pro- 
duced in the short-distance interaction with the same spin, angular mo- 
mentum and colour quantum numbers of the quarkonium state. A single 
parameter, provided by the bound state Bethe-Salpeter wave function or 
its derivative, then describes the hadronization into the observable particle. 
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This model is of course much more predictive: the production rate for a 
colour-singlet QQ state with a given spin and angular momentum, usually 
referred to as QQ[ 2S+1 Lj, 1], can be calculated exactly in pQCD. Moreover, 
the phenomenological parameter can be measured, for instance, in electro- 
magnetic decays and used to make absolute predictions for production rates. 
Despite its physical transparency and predictive power, the CSM is clearly 
not a complete theory: It ignores, at least in its original formulation, rela- 
tivistic corrections which take into account the nonzero relative velocity of 
quark and antiquark in the bound state. Moreover, there are no theorems 
that guarantee that the simple factorization of the quarkonium production 
cross section into a short- distance part and a single non-perturbative param- 
eter is valid in higher orders of perturbation theory. In fact, in the case of 
hadronic production or decay of P-wave quarkonia, the radiative corrections 
to the short- distance cross section contain infrared divergences that cannot 
be factored into a single non-perturbative quantity. This failure is to be at- 
tributed to an incomplete description of the quarkonium wave function, in 
that the Colour Singlet Model ignores the contributions from higher Fock- 
state components. As a result, CSM predictions for producing J/tp and ip' 
states at large pt have been found to grossly underestimate, by more than 
an order of magnitude, the experimental cross sections obtained by the CDF 
collaboration in pp collisions at the Tevatron (see also [B5 and references 
therein). 

The factorization approach by Bodwin, Braaten and Lepage |JT| provides 
a rigorous framework for treating quarkonium production and decays which 
resolves these problems. It extends the CSM by allowing QQ pairs with 
spin, angular momentum and colour quantum numbers different from those 
of the observed quarkonium to hadronize into the latter. In this respect, the 
factorization formalism recovers some of the qualitative features of the CEM. 
The general expression for the production cross section within the NRQCD 
factorization approach reads 

da(H + X)=J2da(QQ[n]+X)(0 H [n]) (1) 

n 

Here da(QQ[n] +X) describes the short-distance production of a QQ pair in 
the colour /spin/angular momentum state n, and (0 H [n\), formally a vacuum 
expectation value of a NRQCD matrix element (see for details), describes 
the hadronization of the QQ pair into the observable quarkonium state H. 
One must note that the cross section is no more given by a single product of a 
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short- distance and a long-distance term as in the CSM, but rather by a sum 
of such terms. Infrared singularities which eventually show up in some of 
the short-distance coefficients would be absorbed into the long-distance part 
of other terms, thereby ensuring a well defined overall result. The relative 
importance of the various contributions in ([!]) can be estimated by using 
NRQCD velocity scaling rules |2C]. For v — > (v being the average velocity 
of the heavy quark in the quarkonium rest frame) each of the NRQCD matrix 
elements scales with a definite power of v and the general expression ([]]) 
can be organized into an expansion in powers of v 2 . For the production of 
S'-wave bound states, like J/ip, the colour-octet matrix elements associated 
with higher Fock state components in the quarkonium wave function are 
suppressed by a factor of v 4 compared to the colour-singlet contribution. 
The NRQCD description of S'-wave quarkonia production or annihilation 
thus reduces to the Colour Singlet Model in the static limit v = O.Q 

However, completely ignoring the octet QQ states, as in the CSM, can 
lead to too small quarkonium production cross sections whenever these states 
can be copiously produced in the short-distance interaction in comparison 
with the colour-singlet ones. This is the case for large-pr J/0 an d produc- 
tion at the Tevatron. Indeed it has been found H that gluon production with 
subsequent fragmentation into a colour-octet 3 Si cc pair (which will even- 
tually hadronize into a physical quarkonium) can successfully describe the 
experimental data, greatly underestimated by the CSM. As will be demon- 
strated below, associated J/ip + 7 photoproduction is a particularly sensitive 
probe of the NRQCD matrix element (O j ^[ 3 Si, 8]) and will therefore pro- 
vide an important consistency check of the colour-octet explanation of the 
high-p T data at the Tevatron. 

A recent review of the factorization approach and its applications to 
quarkonia production at colliders can be also found in ref. 



3 J/i/j + 7 Production 

Associated J/ip + 7 production has been studied in hadronic collisions by 
various authors over the recent years [p3l . Kim and Reya plf pointed out 



that J/tp + 7 photo- and electroproduction represents a powerful tool for dis- 

In the case of P-wave quarkonia, colour-singlet and colour-octet mechanisms con- 
tribute at the same order in v to annihilation rates and production cross sections and 
must therefore both be included for a consistent calculation Ell] . 
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Figure 1: Diagrams contributing within the factorization approach to 
J/ip + photon associated production in direct (a) and resolved (b) photon 
collisions. 

criminating between the Colour Evaporation Model and the Colour Singlet 
Model. The reason for this will be discussed below. Since all previous cal- 
culations were based on the Colour Singlet or Colour Evaporation Model 
it is worth to reconsider J/ip + 7 production within the context of the new 
factorization approach. We will in fact demonstrate that this reaction is a 
particularly sensitive probe of the presence of colour-octet mediated channels 
in J/ip production. 

Before presenting our results within the BBL formalism, let us briefly 
review the reason why associated J/ip + 7 production provides a good dis- 
crimination between the CEM and the CSM in photoproduction. Due to the 
colour constraint of the CSM, the leading order contribution, fig. |TJa, 

7P^7 + S P ^ccf,S 1 ,l]+7^ J/V + 7 (2) 

obviously vanishes in the Colour Singlet Model, since only one gluon is at- 
tached to the cc line. In order to allow for a second gluon one must either go 
to higher order in a s for the direct photon process given above or consider a 
resolved photon process: 

IV 9j + 9 P -^ ccfSi, 1] + 7 -> J/ip + 7 (3) 

In both cases, we can expect a suppression of J/0 + 7 photoproduction within 
the CSM. This is not the case in the CEM. Since the colour constraint is 
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absent, the process can proceed via the direct production channel giving a 
cc pair of invariant mass smaller than that of a DD pair threshold: 

IV -> 7 + 9 P — cc[M c5 < M DS ] + 7 — J/0 + 7 (4) 

Kim and Reya have indeed argued that at HERA within the CEM one should 
expect about a factor of five more J/ip + 7 events than predicted by the CSM. 

As stated before, the factorization approach shares some features of the 
CEM by allowing colour-octet intermediate states. For this case, in par- 
ticular, it allows the J/tp + 7 production to proceed also via direct photon 
interactions through a colour-octet cc pair. To be more specific, the following 
short- distance process can contribute: 

Direct photon (fig. |l|a): 

7 + ft) -^cc[ 1 S ,8]+7 
7 + P ->cc[ 3 S , i,8]+7 
7 + 2 p ^cc[ 3 P/,8]+7 

Resolved photon (fig. |I|b): 

9j + 9 P -^ ccf So, 8] +7 

£7 + 9 P -> cc[ 3 Sa, 1] + 7 (Standard CSM processj 
9i + 9v~^ cc[ 3 Si,8]+7 
9^ + 9p-^ cc[ 3 Pj,8] +7 

Light-quark initiated processes are strongly suppressed and can safely 
be neglected. Resolved production of a colour-singlet xj state and a hard 
photon, with the \Ji radiatively decaying into J/ip plus an unobserved soft 
photon, could constitute a background to our signal, but turns out to be 
zero. 

According to the NRQCD velocity scaling rules PO] and fits to the Teva- 
tron data the colour-octet non-perturbative matrix elements that enter 
associated J/ijj + 7 production should be suppressed by about two orders of 
magnitude compared to the colour- singlet matrix element. Still, however, 
the two following features can be qualitatively expected: 

i) the production of colour-octet states via a direct process - rather than a 
resolved one - will at least partially compensate for the smaller matrix 
elements. We should therefore expect a sizable increase in the overall 
cross section due to colour-octet channels; 
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ii) the distribution in the inelasticity of the J/ip, namely the ratio of its 
energy over the initial photon energy in the proton rest frame, usually 
indicated with z (in a generic frame z = pj/^ ■ p p /p-y ■ p p ), will be more 
peaked towards one for the colour-octet induced processes. This is 
again due to the presence of a direct photon coupling as opposed to the 
resolved one, where the g 1 only carries part of the photon energy into 
the reaction. 

These qualitative features are in fact born out by our numerical analysis to 
be presented and discussed in the next section. The parton cross sections (||) 
have been evaluated using the algebraic computer program FORM [p5[] . We 
note here that the direct photon diagram in fig. |l|a happens to contribute 
to [ 3 5*i,8] production only. Moreover, due to the triple gluon vertex in the 
second diagram of fig. |l]b, the cross sections for the resolved photon contribu- 
tion to 1 Sq and 3 Po,2 octet states display a collinear and infrared singularity 
in the Pt = and z = 1 endpoint. Rather than properly subtracting it we 
will avoid this phase space region by applying a minimum-pr cut. 



4 Numerical Results 

The numerical results have been obtained by adopting the following set of 
parameters. As for the strong coupling, we employ the leading order formula 

asifi) = (33 - 2^k^m (5) 

with n/ = 4 and A£q = 120 MeV. The scale fx has been taken equal to the 



transverse J/ip mass, namely fx = M T = \j^rn^ c + pj,, and the charm mass 
m c has been taken equal to 1.5 GeV. For both the proton and the resolved 



photon structure functions the GRV92 leading order parton densities p6 
have been employed, evaluated at the factorization scale fxp = M T . 

Finally, the non-perturbative matrix elements have been fixed to values 
compatible with both the NRQCD scaling rules and with the Tevatron fits, 
as in ref . M . Their values are summarized in Table [I]. 

We first compare the total cross sections given by the various production 
channels, with a fixed px > 1 GeV cut applied to avoid the infrared/collinear 
singularities in some of the cross sections. Table shows the results for a) 
7P collisions in a fixed target set-up, with an incoming photon energy of 100 
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{0 J W[ 6 S u l\) 


1.16 GeV 3 




(0 J ^[%,8]) 
(O j ^[ 3 S u 8}) 
(O J ^[ 3 P ,8])/m 2 c 


10- 2 GeV 3 


1 7 


10- 2 GeV 3 


S 7 


10" 2 GeV 3 


S 7 



Table 1: Values of the NRQCD matrix elements used in the numerical 
analysis, with the velocity and mass scaling, v is the velocity of the 
heavy quark in the quarkonium rest frame. For charmonium it holds 
v 2 ~ 0.23. It also holds {O j H> [ 3 Pj,8]} ~ (2J + 1)(O j ^ [ 3 P ,8]). 



GeV; b) jp collisions at a center of mass energy of 100 GeV and c) HERA 
kinematics, i.e. 27.5 GeV electron and 820 GeV proton collisions. For the 
case c) the electroproduction cross sections have been evaluated by averaging 
the photoproduction ones weighted by the usual Weizsacker- Williams flux 
factor, 
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(7) 



E 7 /E e , Qj 



m 2 e y/(l — y) and m e is the electron mass. We 



4 GeV 2 and y m in = 0.15, y ma x = 0.86 according to [^7 



where y = 
adopt Q 2 m 

The results presented in Table ^| clearly show that in a typical fixed target 
set-up, with a photon beam energy of 100 GeV, the cross section is domi- 
nated by the direct photon production of a colour-octet 3 Si state. Although 
the cross section is small, the observation of this process at fixed-target ex- 
periments would therefore already provide good evidence for the importance 
of colour-octet contributions. 

At a higher center of mass energy, on the other hand, the production of 
a Jfif) via a colour-singlet 3 Si state in resolved photon collisions and via a 
colour-octet 3 Si state in direct photon collisions represents the major part of 
the cross section. The numbers presented in Table |2] indicate that at HERA 
energies colour-octet channels could amount to about 40% of the overall 
production rate. The presence of colour-octet contributions can however not 
be assessed from total cross sections alone, given the large normalization 
uncertainties present in the calculation from higher order corrections, parton 



8 







IP 


IP 


ep 


ep 


Channel 


-fc 7 = 100 GeV 
PT > 1 GeV 


y/s = 100 GeV 
Pt > 1 GeV 


= 300 GeV 
Pt > 1 GeV 


sj s = dOO GeV 
p T > 1 GeV 












-3 < 77 < 3 












E-f > 2 GeV 




•S'o, 8 










Direct 




.48 


7.67 


.63 


.44 




3 Pj,S 












So, 8 


.0013 


.35 


.044 


.020 




s 5i,l 


.072 


16.70 


2.05 


.75 




3 ^i,8 


.0012 


.27 


.033 


.012 


Resolved 


3 Po,8 


.0046 


1.03 


.13 


.048 




3 Pi,S 


.0005 


.14 


.018 


.007 




3 P 2 ,8 


.0045 


.97 


.12 


.043 



Table 2: Results for the total cross sections (in pb). Notice that while 
in the first data column £L refers to the incoming photon beam energy 
in the last one it refers instead to the produced outgoing photon. 



distribution functions and charm quark mass values, as well as unknown 
higher-twist contributions. 

Therefore, we propose to study differential distributions to disentangle 
the colour-octet contributions from the standard colour-singlet one. 

In fig. |2] we show the differential distributions related to the total jp cross 
sections at y/s = 100 GeV with a minimum-p T cut of 1 GeV, presented in 
Table |2|. The distributions due to colour-singlet [ 3 Si,I] production in resolved 
photon collisions (continuous line) and to colour-octet [ 3 Si,8] production in 
direct photon collisions (dashed line) only are shown. The distributions due 
to the other colour-octet processes do indeed present the same features as 
the ones of [ 3 Si,l], being also produced in resolved photon interactions, but 
are suppressed in magnitude, as can be seen from Table |2|. Their inclusion 
would therefore not change the picture we are going to discuss. 

As expected, the effect of the colour-octet [ 3 «Si,8] contribution produced in 
direct photon processes can easily be seen in at least some of the plots. While 
the pt of the J/ip and the invariant mass distribution M of the J/0-7 pair 
are pretty similar for the colour- singlet and colour-octet induced channels, 
the z, rapidity and photon energy distributions do indeed show a strikingly 
different behaviour. 
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Figure 2: Differential distributions in jp collision at y/s = 100 GeV. A 
minimum-pr cut of 1 GeV is applied. 



10 



Recalling that we put ourselves in the so-called "HERA- frame" , with the 
photon (or the electron) traveling in the direction of negative rapidities, we 
notice how the direct photon coupling favours the production of the quarko- 
nium and of the photon in the negative rapidities region. This contrasts the 
case of resolved photon production of colour-singlet 3 Si states, which are 
uniformly produced around the central rapidity region. 

As for the z distribution, the resolved photon process predicts a decrease 
of the cross section going towards the high-z region.. The direct photon 
process does on the other hand predict the opposite behaviour: the cross 
section now increases going towards z — 1. The small dip in the last few bins 
is due to the minimum-p^ cut. 

The photon energy distribution behaves similarly to the z distribution, 
and is predicted to be much harder in direct photon processes. 

The above distributions (the shapes of which have been checked to be 
robust with respect to a higher px cut, to be sure of the absence of p™ m 
effects) provide a clear experimental signature; in particular the observation 
of a substantial fraction of J/ip + 7 events in the high-z region would already 
be good evidence for the presence of colour-octet contributions to the J/i/j 
production cross section. 

However, since actually HERA provides electron-proton interactions, we 
are going to investigate how the distributions look like in this frame. For 
consistency, we will first consider, in fig. |3|, the differential cross section cor- 
responding to the total ep cross section shown in Table Q Successively, the 
effect of some experimental-like cuts will be taken into account. 

Comparing fig. |] with fig. || we notice that, due to the system being now 
boosted in the proton direction (i.e. positive rapidities), the non Lorentz 
invariant observables are affected. More precisely, the rapidity distributions 
of both the J/ip and of the photon are slightly smeared and shifted towards 
positive rapidities. The largest difference can however be observed in the 
photon energy distributions: the one due to the resolved photon process is 
hardened with respect to 7p interactions, while the one due to direct photon 
interaction is greatly softened. This can be understood assuming that most 
of the very energetic photons in direct photon collisions are produced in the 
incoming photon direction (see the rapidity distribution). The effect of a 
boost of the event in the opposite direction will then soften this distribution. 
In resolved photon interactions, on the other hand, events of this kind are 
produced more uniformly around the central rapidity region, and the boost 
will then also harden at least part of them. 
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Figure 3: Differential distributions in ep collision at ^fs = 300 GeV. A 
minimum-pr cut of 1 GeV is applied. 
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Figure 4: Differential distributions in ep collision at v / s = 300 GeV. A 
minimum-j>r cu t of 1 GeV, a minimum outgoing photon energy cut of 
2 GeV and a pseudorapidity cut \rjj/^^\ < 3 are applied. 
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The effect of applying some realistic experimental cuts can be appreciated 
in fig. f|. Since the produced photon must be clearly visible in the event, it is 
now required to lie (together with the J/ip) within the pseudorapidity region 
\t)\ < 3 and to have an energy, in the lab frame, greater than 2 GeV. This last 
requirement also ensures that it will not be mistaken with a photon coming 
from the radiative decay of a P-wave state to a J/ip, having energies of a few 
hundreds MeV. A further experimental selection criterion to exclude radiative 
decays photons consists in asking the photon px to be roughly opposite to 
that of the J/ip. The distributions of fig. (| are qualitatively similar to the ones 
- without cuts - of fig. |3], showing that colour-octet contributions in associated 
J/0+7 electroproduction should be visible in the photon fragmentation region 
(negative rapidities in the HERA frame) or, more clearly, in the large- z 
region. 

5 J/0 + 7 in hadron collisions 

Before closing, it is worth noticing that the resolved processes of fig. |I]b also 
contribute to the associated J/ip + 7 production in hadron collisions. To 
check whether a signature for colour-octet mediated processes can be seen 
at a hadron collider we have evaluated the total J/ip + 7 cross section in pp 
collisions at y/s = 1800 GeV. Light-quark initiated processes are suppressed 
and have been neglected. The following cuts have been imposed: 

p T > 4 GeV, Irjj/^l < 0.6, E 7 > 2 GeV (8) 

The results are displayed in Table |||. The colour-octet terms can be seen to 
enhance the cross section by about 50%. However, given the large normaliza- 
tion uncertainties involved and the fact that the differential distributions in 
this case do not differ substantially from the colour-singlet process' ones, we 
conclude that there is little hope to shed light on the colour-octet mechanism 
via J/ip + 7 production in hadron collisions. 

6 Conclusions 

The associated J/0 + 7 production in ep and 7p collisions has been proposed 
as a powerful tool for establishing the presence of quarkonia production pro- 
cesses mediated by colour-octet QQ states, as suggested by the factorization 
approach of Bodwin, Braaten and Lepage. 
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Vpp{J/ip + V (pb) 


Channel 


y/s = 1800 GeV 




Pt > 4 Lrev 




kj/^l < °-6 




^ > 2 GeV 


*S u l 


93.2 


So, 8 


12.3 


3 S U 8 


1.5 


3 Po,8 


10.0 


3 Pi,8 


13.8 


3 P 2 ,8 


7.7 



Table 3: Results for the total cross sections in hadron collisions (in 
pb). 



The very fact that colour-octet contributions to J/0 + 7 photoproduction 
can proceed also via a direct photon coupling rather than via a resolved 
one only - as for the colour-singlet channel - leads to clean experimental 
signatures. We have shown that at HERA energies the colour-octet terms 
can increase the cross section by about 50% and, most importantly, produce 
a J/ip energy distribution da/dz strikingly different from the one predicted 
by the colour- singlet channel alone. 

The hadroproduction case, taking the Tevatron as an example, has also 
been investigated. This cross section is also increased by about 50% by octet 
terms, but no significant signatures in differential distributions can be found. 



Note added: After completion of our work an analysis of colour-octet 
contributions to associated J/ip + 7 production in hadron-hadron collisions 
appeared [28]| . After correcting what is probably a trivial misprint in the 
labelling of figs. 6-8 of |28|] (they should have units nb rather than pb), their 
results are consistent with ours. 
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